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ABSTRACT 


The  first  part  of  this  paper  presents  an  algebraic 
topological  approach  to  finding  all  paths  in  a  given 
graph  under  various  conditions.  If  a  path  with  any 
particular  characteristics  is  desired,  such  as:  mini¬ 
mum  distance,  least  cost,  most  reliable,  etc.,  it  can 
be  extracted  from  the  list  of  all  paths.  In  the  second 
part  of  the  paper,  a  direct  method  for  finding  a  mini¬ 
mum  path  is  given.  Both  techniques  are  straight-forward 
and  can  easily  be  performed  by  a  computer. 
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INTRODUCTION 

This  paper  treats  two  problems  related  to  graphs  (i.e.  diagrams 
consisting  of  nodes  and  branches).  The  first  problem  is  to  find  all 
of  the  paths  in  a  graph,  from  any  one  node  to  another.  The  second 
problem  is  to  find  the  minimum  length  path  In  a  graph,  from  one  node 
to  another.  The  mathematical  techniques  Involved  are  matrix  algebra 
and  topology. 

In  the  first  part  of  the  paper  the  general  path  finding  problem 
is  handled.  At  first  the  theory  Is  presented  and  then  detailed  ex¬ 
amples  are  given.  In  the  second  part  of  the  paper  the  minimum  path 
problem  Is  handled.  Similarly,  the  theory  Is  presented  at  first  and 
then  detailed  examples  are  given.  Detailed  proofs  of  all  theorems 
are  to  be  found  in  the  appendix. 
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PART  I  -  GENERAL  PATH  FINDING 
Topological  Foundations: 

This  part  presents  a  way  of  finding  all  paths  from  one  point  in  a 
graph  to  any  other  point.  This  technique  is  the  elaboration  of  ideas 
presented  in  O-tasc*.  Various  restrictions  can  be  placed  upon  the  paths 
so  that  they  include  certain  branches  and  do  not  include  certain  other 
branches . 


Let  us  assume  a  graph  to  be  given  and  hence  its  node-branch  inci¬ 


dence  table.  As  an  example: 


1 

2 

3 

4 

5 

6 

A 

1 

1 

0 

0 

0 

1 

B 

1 

0 
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0 

C 

0 

1 

1 

0 

1 

0 

D 

0 

0 

0 

1 

1 

1 

Incidence  Table 


Fig-  1 


where  a  one  in  the  ij  positions  means  that  node  i  is  incident  to  branch 
j  and  a  zero  means  that  it  is  not  incident.  For  the  details  of  such 
topological  concepts  see  V-AS.  The  following  is  a  brief  summary  of  the 
1 -dimensional  incidence  table  and  certain  of  its  properties.  Let  us  de¬ 
note  a  branch  or  set  of  branches  by  its  component  vector  expression. 

Thus  in  the  above  figure  the  set  of  branches  1  and  3  would  be  denoted  by 
(101000),  where  the  l's  in  the  first  and  third  positions  mean  that 
branches  1  and  3  are  present  and  the  0's  in  the  second,  fourth,  fifth 
and  sixth  positions  mean  that  branches  2, 4, 5, 6  are  not  present.  Similarly 
the  set  of  nodes  A  and  D  are  denoted  by  (1001).  The  boundary  of  a 
branch  is  defined  to  be  the  nodes  to  which  this  branch  is  incident.  The 
boundary  of  branch  1  is  written 


B  (100000)  -  (1100) 


(1) 


*  The  above  abbreviations  or  acronyms  have  been  found  to  be  a  very  useful 
way  of  referencing  the  literature.  The  acronyms  are  listed  in  alphabeti¬ 
cal  order  in  the  reference  section,  with  the  full  reference  attached. 

The  part  of  the  acronym  before  the  hyphen  represents  the  author  and  the 
part  after  the  hyphen  the  title.  Capital  letters  following  the  hyphen 
signifies  a  book,  and  small  letters  a  journal  article;  if  only  the  first 
letter  is  capitalized,  then  this  signifies  an  individually  available  report. 
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where  S  is  the  boundary  operator.  The  co-boundary  of  a  node  is  defined 
to  be  the  branches  which  are  incident  to  this  node.  The  co-boundary  of 
node  A  is  written: 

6  (1000)  -  (110001)  (2) 

where  5  is  the  co -boundary  operator.  Thus  the  rows  and  columns  of  the 
incidence  table  correspond  to  the  coboundaries  of  the  nodes  and  boundaries 
of  the  branches. 

Let  us  define  an  additive  operation  between  two  binary  component  vec- 


.  .a  )+(b 
n  I 

•bn>  ■ 

(Cj,... 

c  )  in  the 
n 

following  manner: 

i) 

a.-o, 

b^=o 

then 

c  .  =o 
i 

u) 

V1, 

o 

ii 

43 

then 

c .  =1 
i 

iii) 

a .  =o , 
l 

b  *1 
i 

then 

V1 

iiii) 

V1* 

b  «1 
i 

then 

c^*o 

this  can  be  considered  as  a  modulo  2  position-wise  addition. 

Let  us  define  the  boundary  of  a  set  of  branches  to  be  the  sum  of  the 
boundaries  of  each  branch  in  the  set.  Let  us  define  the  coboundary  of  a  set 

of  nodes  to  be  the  sum  of  the  coboundaries  of  each  node  in  the  set. 

Let  us  define  a  scalar  multiplication  of  vectors  by  the  scalars  o  and  1 

as  zero  times  a  vector  is  the  zero  vector  and  1  times  a  vector  is  the  vector 

itself . 


It  is  then  quite  natural  to  view  the  incidence  table  as  a  matrix  with 
addition  and  multiplication  between  rows  and  columns  as  defined  above.  It 
can  be  shown  (see  V-AS)  that  these  rows  form  a  vector  space  under  the  field 
of  integers  modulo  2,  whose  rank  is  equal  to  the  number  of  rows  minus  1.  It 
can  also  be  shown  that  the  boundary  and  coboundary  operations  can  be  re¬ 
placed  by  matrix  multiplication  involving  the  incidence  matrix,  as  exempli¬ 
fied  in  equations  3  and  4, 


(3) 


(4) 


Let  us  define  two  branches  ^^>^2  to  c°lncl^ent  if  b^  and  b^  are  in 
cident  tc  the  same  node  at  exactly  one  of  their  ends. 
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0 

= 
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_ 
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1 

noooi 

r  n 

r 

1 

101100 

r 

1 

6(1000)» 

1000 

H, 

* 
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• 
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• 

110001 

L  J 
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000111 

L 
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Let  us  define  a  set  of  branches  b^,...bn  to  be  a  path  from  a  node  A 
to  a  node  B,  if  bj  is  incident  to  node  A  and  is  coincident  with  bj^ 
(i=l , 2. . ♦ n-1) ,  and  bn  is  incident  to  node  B.  Any  permutation  of  the 
set  of  branches  in  a  path  will  also  be  considered  as  a  path. 

Theorem  1 :  The  boundary  of  a  path  b^,...b^  from  node  A  to  node  B  is  the 
set  of  nodes  A  and  B. 

Proof:  See  Appendix. 

Path  Calculation 


If  a  graph  is  given  (i.e.  its  incidence  matrix  H)  and  there  exists  a 
path  p  from  node  A  to  node  B  (let  us  denote  the  vector  expression  of  these 
two  nodes  as  n)  then 


H  •  p  =  n  (5) 

Also,  any  path  x  from  node  A  to  node  B  must  satisfy 

H  •  x  =  n  (6) 

Thus  the  general  problem  of  finding  a  path  p  from  node  A  to  node  B  is  to 
solve  the  above  matrix  equation  for  x. 

The  following  theory  presents  a  direct  method  for  finding  the  general  so¬ 
lution  to  equation  6.  If  the  number  of  rows  in  H  is  a  then  it  is  easily 
shown  (see  V-AS)  that  the  rank  of  H  is  Ol  and  any  one  of  the  equations  ob¬ 
tained  from  (6)  is  linearly  dependent  on  the  rest  and  can  be  eliminated.  Let 
us  call  H'  the  matrix  obtained  by  deleting  one  of  the  rows  H  and  n'  as 
the  column  matrix  obtained  by  deleting  the  corresponding  entry  in  n.  Then 
(6)  is  equivalent  to: 


H1-  x  =  n'  (7) 

Definition:  A  rectangular  matrix  of  size  r  x  s,  where  s  >  r,  will  be 
called  semi -diagonalized  if  its  left  r  x  r  part  is  a  unit  matrix  and  its 
right  r  x  (s-r)  part  is  an  arbitrary  matrix  of  zero's  and  one's  denoted  by 
H  . 

c  —I 


c 


1 


r 


s-r 
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Theorem  2:  The  matrix  H1  can  be  brought  into  semi -diagonalized  form 
by  using  the  following  two  operations: 

I  The  addition  of  one  row  to  another 

II  The  interchange  of  two  columns 

Proof :  See  Appendix 

Theorem  3:  If  H*  is  the  semi-diagonalized  form  of  H'  and  x*  is  ob¬ 

tained  from  x  by  the  corresponding  operations  of  (II)  and  n*  is  obtained 
from  n'  by  the  corresponding  operations  of  (I)  then  the  solutions  to  the 
matrix  equation 


H* .  x*  «  n* 


(3) 


are  equivalent  to  equation  (7), 

Proof:  See  Appendix 

The  rest  of  the  calculation  will  be  developed  by  the  use  of  matrices 
directly. 


Let: 


H*  = 


(9) 


„  *  * 
x*  =  HXj^  .  .  .  xg 


n*  =  |jn1*  .  .  .  nr* 


(The  matrix  nctatio 
From  Equation  (8) 

r 


n[ 


means  transpose) 


“  1 


H  * 
c 


X  * 
1 


X  * 

a 


V 


n  * 
r 


(10) 


(11) 
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which  is  equivalent  to 


r—  — ■* 

—  — 

_  — 

r 

—  — 

1 . 

x  * 
X1 

+ 

H  * 

x  * 
r+1 

“i* 

1 

x  * 

r 

c 

x  * 

n  * 

—  — 

_  s  _ 

! 

r  _ 

It  then  follows  that 


-  — 

r_  - 

— • 

—  — 

x  * 
X1 

ni* 

+ 

H  * 

X  * 

r+1 

X  * 

n  * 

c 

x  * 

r 

r  _ 

—  — 

_  s  _ 

(12) 


(13) 


due  tc  the  modulo  2  properties  of  these  matrices.  If  (xr+l*>  •••  xs*) 
is  allowed  to  vary  through  all  binary  (s-r) -tuples  from  (0,0,... 0)  to 
(1,1 ,...1)  then  all  solutions  for  x  are  obtained*  The  solution  can  be 
brought  into  a  more  elegant  form  by  enlarging  both  sides  in  the  following 
manner : 


x  * 

X1 

x  * 

X 

H  * 
c 

X  * 

r+1 

r 

x  ,, 
r+1 

x  * 

s  _ 

r 

‘o' 

_0  _ 

+ 

1 

1  _ 

x  * 
s 

(14) 

If  both  sides  are  re-ordered  by  the  appropriate  interchange  of  rows  and 
columns  on  both  sides  of  the  equation  so  as  to  obtain  x  in  its  original 
order  on  the  left-hand  side,  then 


x  =  y  +  H  •  x 1  (15) 

c 

where  y  is  a  particular  solution  for  x,  and  depends  only  upon  the  n 
columns  of  H*  which  were  used  in  the  semi-diagonalization  (which  corres¬ 
pond  to  a  tree  of  the  graph);  Hc  depends  only  on  the  semi-diagonalized 
Hc ' ;  x1  is  a  subset  of  the  x  and  corresponds  to  the  choice  of  co-tree 
(complement  of  a  tree,  or  the  s-r  columns  not  in  the  semi-diagonalized 
part  of  H  *). 
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The  total  number  of  solutions  for  x  can  easily  be  obtained.  These 
solutions  were  found  by  allowing  (xr+]** • *x8*)  to  vary  through  all  s-r 
tuples  of  binary  numbers.  Hence  there  are  2(8“r)  solutions  for  x.  If 
we  let  n  be  the  number  of  nodes  in  the  graph ,  b  be  the  number  of  branches 
in  the  graph,  then  r  ■  n-1  and  s  -  b.  Hence  the  number  of  solutions  K 
is  : 

K  -  2(b'n+1^  (16) 


Excessive  Solutions: 

In  the  previous  section  we  showed  how  to  solve  the  system  of  equa¬ 
tions  in  (6).  In  addition,  we  showed  that  any  desired  path  will  satisfy 
(6)  and  hence  be  included  in  the  general  solutions.  It  has  not  been  shown 
that  any  element  of  the  general  solution  will  be  a  desired  path.  This 
section  will  serve  to  clarify  this  last  point. 

Let  us  define  a  loop  as  a  set  of  three  or  more  branches  b  ,...b  such 
that  b^  is  coincident  to  b  , i=l , 2 , . , . n  and  b^  is  coincident  to  b^, 
where  no  subset  of  the  b.  satisfy  the  same  condition.  A  set  of  two  branches 
both  of  whose  ends  are  tfte  same  will  also  be  called  a  loop.  Any  branches 

both  of  whose  ends  are  the  same  will  also  be  called  a  loop.  Any  permuta¬ 
tion  of  a  set  of  branches  which  form  a  loop  will  be  considered  a  loop. 

Theorem  4 :  If  b  is  a  set  of  branches  then  the  equation 

H  •  b  =  0  (17) 

is  satisfied  if  and  only  if  b  is  a  loop  or  a  set  of  loops. 

Proof :  See  Appendix 

Definition :  If  two  vectors  of  zero's  and  one's 

a  =*  (a^, . .  .  a^) 
b  =  (bl,...bR) 

are  given,  where  a  £  b  and  if  a.  =  1  implies  that  b. 
be  said  to  be  included  in  b,  written  a  ch.  1 

The  following  theorem  completely  cla  ^ifies  the  solutions  to  equation 

(6). 

Theorem  5 :  If  b  is  a  set  of  branches  and  n  is  a  set  of  two  distinct  nodes 
A  and  B ,  then 

H  .  b  =  n  (19) 

if  and  only  if  either  b  is  a  path  from  A  to  B  or  b  is  the  sum  of  a  path 
k  from  A  to  B  and  loops  , . , . where  k  b ,  i^czb,  ...  ?^cb. 

Proof:  See  Appendix 


(10) 

=  1,  then  a  will 
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Thus  the  solutions  to  equation  (6)  are  of  two  kinds;  the  paths 
which  were  required,  and  the  paths  plus  excessive  loops.  Once  all  the 
solutions  have  been  found,  there  are  several  ways  of  eliminating  those 
solutions  which  contain  excessive  loops  (excessive  solutions). 

The  first  method  is  based  upon  the  partial  ordering  induced  upon 
the  solutions  by  the  above  defined  inclusion  relation. 

Theorem  6 :  If  a  set  of  branches  b  which  is  a  solution  to  equation  (6) 
contains  loops  ^  cb,...i  cb  then  there  exists  a  path  p 

such  that  p  is  also  a  solution  of  equation  (10)  and  p  cb. 

Proof :  See  Appendix 

Thus  the  method  for  eliminating  excessive  solutions  is  as  follows: 
find  all  of  the  solutions,  then  all  solutions  b  such  that  there  exist 
smaller  solutions  p  (i.e.  pcb)  are  excessive  solutions  and  should  be 
eliminated. 

A  second  method  to  eliminate  excessive  solutions  is  by  rank  considera¬ 
tions  and  can  be  found  in  O-tasc  276-9. 

Routes  Under  Specified  Conditions: 

In  general  it  is  possible  to  specify  that  certain  branches  should 
be  included  in  the  path  and  that  certain  other  branches  should  not  be  in¬ 
cluded,  In  the  original  formulation  H«x  =*  n  the  components  of  x  were 
all  assumed  to  be  unknown.  Let  us  now  assume  that  x  will  be  partially 
specified.  Simply  set  equal  to  zero  those  branches  which  should  not  be 
included  and  set  equal  to  one  those  branches  which  should  be  included.  If 
it  is  desirable,  the  general  solution  of  equation  (6)  can  first  be  found 
and  then  these  conditions  can  be  used  as  constraints. 

Examples : 

Some  illustrative  examples  will  help  to  clarify  many  of  these  ideas. 
Figure  2  is  an  imaginary  telephone  communication  network.  The  matrices  in 
Figure  3  are  then  obtained. 


10,026 


I  2  3  4  5  6  7  6  9  10  II  12  13  14 

CHICAGO  A  I  I  I  I  I 

BOSTON  B  I  II 

WASHINGTON  Cl  II 

ALASKA  D  ’ 

SAN  FRAN.  E 

COL.  SP.  I  F 

GREEN.  G 

NEW  YORK  H 

C.C.  I 

LOS  ANGLES  J 

C0L.SP.2  K 


INCIDENCE 

MATRIX 


I  2  3  4  5  6  7  8  9  10  II  12  13  14 


I  I  I 

I  I  I 

I  I  I 

I  I  I 

I 

I  I 

I  I 

I  I 

I  I 
I  I 


REDUCED 

INCIDENCE 

MATRIX 

H’ 

(A-REMOVED) 


B+G+H 

C+I 

D 

E+J 

F 

G 

H 

l+K 

J 

K 


I  2  3  4  5  6  7  9  12  13  B  10  II  14 


I  I  I 

I  I  I 

I  I  I 

I  I  I 


SEMI -DIAGONALIZED 
MATRIX 


TM-  3421 


11 


Problem  1 

Let  us  find  all  paths  from  Alaska  to  New  York.  This  is  the  same  as 
finding  the  general  solution  to 


it  then  follows  that 


where 


H.x  -  |l 
t 

H.x*  -  £ 


00010001000 


1010001000 


(23) 


(24) 


and 


x*  =  (x1)x2,x3>x4,x5)x6)x7>x9>x12jx13.x8>x10>x11)x14) 


(25) 
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C11 

C14 


(26) 


The  solutions  are  easily  found  by  substituting  all  binary  values  for 

(VX10’X11>X14): 
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The  solutions  marked  with  an  x  mean  those  with  excessive  loops.  The 
sixteen  diagrams  on  pages  19  to  22  give  the  actual  sixteen  paths. 


Problem  2 

In  this  problem  let  us  assume  that  in  Figure  1  Chicago  has  been  des 
troyed  and  that  we  wish  to  find  all  paths  from  Alaska  to  New  York.  One 
way  to  find  the  solution  would  be  to  start  from  the  beginning  and  define 
a  new  incidence  matrix  with  Chicago  removed.  However,  a  much  shorter  me 
thod  is  possible  starting  from  the  solution  to  Problem  1  in  Equation  26. 
The  above  condition  can  be  specified  by  forcing: 


(x1,x2,x3,x^,x5)  =  (0,0, 0,0,0) 
in  Equation  25.  This  gives 


(23) 


(29) 
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an  independent  sub -set  of  these  equations  is 


1 

1010 

x8 

0 

1001 

• 

*10 

1 

0110 

X11 

—  — 

—  — 

x14 

and  the  solution  is  easily  found  as 


X. 

0~ 

r 

8 

x10 

0 

i 

X11 

1 

+ 

i 

X14 

0 

i 

—  _ J 

— 

— 

(30) 


(31) 


If  this  is  substituted  in  Equation  26  we  obtain  the  following  two  solu¬ 
tions  : 


i. 

pi 

=  Wxn 

2. 

p2 

=  x8,x9,x10 

(32) 


Problem  3 

Let  us  find  all  paths  from  Alaska  to  New  York  which  pass  and  do 
not  pass  x  .  Again  we  can  start  from  Equation  26.  From  Equation  26  we 
obtain 


0 

_ 

1 

+ 

0110 

1 

0 

1000 

—  —J 

_  _ 

The  solution  is  easily  found 


"10 

hi 

C14 


(33) 


x8 

100 

1 

xio 

no 

X11 

X11 

010 

X14 

X14J 

001J 

(34) 


If  this  is  substituted  in  Equation  26  we  obtain  the  following  four  solu¬ 
tions  : 
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(35) 


Part:  II  -  Determining  Minimum  Paths 

In  the  study  of  networks  the  problem  of  determining  the  shortest  path  (in 
the  network)  joining  two  given  points,  arises.  This  problem  may  be  generalized 
to  determining  the  least  cost  paths  connecting  two  given  points.  If  a  non¬ 
negative  number  s  assigned  to  each  branch  of  a  network,  then  to  each  path  in 
the  network  a  number  called  its  cost,  may  be  assigned.  This  number  is  the  sum 
of  the  numbers  assigned  to  the  branches  comprising  the  path.  Various  algorithms 
exist  for  finding  the  least  cost  paths  connecting  two  points  of  a  network 
(p3-T0G],  Ch  .7,  (TT-apaj  M-SPTM  ,Vol.  2).  We  shall  present  a  matrix  method  which 
is  applicable  to  both  directed  and  undirected  networks.  In  the  main  body  of 
this  paper  we  shall  formulate  the  method  for  undirected  networks,  and  in  an  ap¬ 
pendix,  the  method  will  be  formulated  in  general. 


Let  N  be  a  network  with  a  finite  number  of  branches  and  vertices  (  nodes). 

We  say  that  two  points  of  N  are  joined  by  a  branch  if  these  points  are  the  end 
points  of  the  branch.  We  shall  assume  that  no  point  of  N  is  joined  to  itself 
by  a  branch.  By  a  cost  function  on  N  we  mean  a  function  f  that  assigns  to  each 
pair  of  vertices  in  N,  either  a  non-negative  real  number  or  the  indeterminate 
quantity  «>.  Furthermore  f  obeys  the  conditions 

(1)  f(v,v)  =  o  for  all  vertices  v 

(2)  f(v,w)  =  f(w,v)  for  all  vertices  v,w 


(This  condition  will  be  dropped  later). 

(3)  If  v^w  then  f(v,w)is  a  non-negative  real  number  if  v,w  are 
joined  by  a  branch  of  N, 

(4)  If  v  w  then  f(v,w)  =  »  if  there  is  no  branch  joining  v  to  w. 


The  indeterminate  quantity  »  obeys  the  conditions: 

(1)  oo  +  r  =  oo  for  any  real  number  r. 

(2)  00  +  00  =  00 

(3)  oo  >  r  for  any  real  number  r. 

Let  the  vertices  of  N  be  v1  ,v ^  .  .  .Vr.  Then  with  N  we  associate  the 
symmetric  matrix 
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is 


f(vl,vl>  •  •  •  •  f (vj  >VR) 


f(v  V,),.  .  .  .  f(v  ,v  ) 
n  i  n  n 


Let  us  consider  an  example  of  a  network  and  a  cost  function  on  that  network. 


The  matrix  associated  with  this  network  and  cost  function  is: 


v 

v 

v 

v 

v 

v 


1 

2 

3 

4 

5 

6 


00 


5 

00 

1 


1  «  5 
0  5  oo 
5  0  2 
oo  2  0 
oo  4  5 
00  00  2 


00 


1 


4  00 

5  2 

0  9 

9  0 


This  matrix  not  only  tells  us  what  cost  is  assigned  to  each  branch,  but  in 
addition  tells  us  what  vertices  are  joined  by  branches.  For  instance,  we 
know  and  are  not  joined  by  a  branch  since  the  (1,5)  entry  is  ». 
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This  follows  from  conditions  2  and  3  for  cost  functions. 

We  shall  now  define  a  binary  composition  on  matrices  whose  entriea  are 
real  numbers  or  the  indeterminate 

Definition: 


— 

— 

— 

—  — 

8ir  • 

*  4 *ln 

* 

V  • 

■bln 

m 

c. , .  .  .c, 

11  In 

8nl  *  ' 

.  .  a 

nn 

bnl-  ‘ 

*  b 

nn 

i 

c  . .  .  .c 
nl  nn 

— 

__ 

— 

_ 

_  _ 

where  =  Min  («u  +  »12  +  »>2j . »in  +  ) 

Definition:  Let  A  be  a  matrix  whose  entries  are  either  real  numbers 
or  the  indeterminate  ».  Then  A'*'  *  A,  A^2'  -  A  *  A,  A^'  -  A^  '  *  A  .  .  . 

An  -  A^11-1)  *  A. 

Lemma  1:  Given  matrices  A,B,C 

A*(B*C)  -  (A*B ) *C 

Lemma  2:  Given  matrix  A, 

A0 0  *  A(k')  w  A(k+k') 

Again  let  A  be  the  matrix  associated  with  a  given  network  N  and  let 
f  be  a  cost  function  on  N.  We  shall  see  that  the  entries  of  A^n*^,  where  n 
is  the  number  of  vertices,  are  the  minimum  costs  of  getting  from  one  vertex 
to  another.  Consider  the  matrix  A*A  -  a(2).  Let  a^'  be  the  i,jth  entry 
in  A<2>.  Then  *£j>  is  Min  (  *n+*x  j  >ai2+A2  j  *  ain+*nj^Where  aij"f *vi,vj'- 

The  cost  of  a  least?cost  path  between  v  ,v  consisting  of  no  more  than  2  branches 
is  thus  given  by  a^j  .  Similarly  the  following  theorem  can  be  shown. 

Theorem  7.  Let  a^^  be  the  i,jth  entry  of  A^.  Then  a^  is  the  least  cost 
entailed  by  any  path  of  <  k  branches,  joining  v^  to  v j  .  J 

Now, since  a  least  cost  path  has  no  loop  in  it,  it  suffices  to  consider  the 
matrix  A^n"^  to  determine  the  least  cost  to  get  from  one  point  to  another. 
Clearly,  any  path  having  n  branches  or  more,  must  have  a  loop  so  our  assertion 
follows  from  Theorem  7,  and  condition  3  for  a  cost  function. 

(k) 


Theorem  8:  If  k'<k,  then  each  entry  of  A 
corresponding  entry  of  A(k’). 


is  greater  than  or  equal  to  the 
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Thus  for  k>n-l,  A^-A^n_1\ 


Let  us  now  consider  the  network  and 
coat  function  we  gave  as  an  example  (See  Fig.  1).  Let  A  be  their  asso¬ 
ciated  matrix.  A  is  a  6x6  matrix  so  it  suffices  to  consicj*F  A'*', 

XAW* 

so 

we  need  3  compositions  instead  of  4.  Here  are  the  compositions  done  in 
detail : 


ciet&a  matrix,  a  is  a  o  x  o  matrix  so  it  surricea  to  consider 

But  Lemma  2  gives  us  a  simpler  way  of  getting  A^.  Notice 

A(4),  A(4)  a  a(2'*A'2',  and  AW^A*A.  By  the  above  remark  A  (3)*A^' , 


(2) 


0 

1 

00 

5 

00 

1 

0 

1 

7 

3 

10 

1 

1 

0 

5 

00 

00 

00 

1 

0 

5 

6 

9 

2 

00 

5 

0 

2 

4 

00 

« 

6 

5 

0 

2 

4 

4 

5 

00 

2 

0 

5 

2 

3 

6 

2 

0 

5 

2 

00 

00 

4 

5 

0 

9 

10 

9 

4 

5 

0 

7 

1 

00 

00 

2 

9 

0 

1 

2 

4 

2 

7 

0 

(2) 

0 

1 

6 

3 

10 

1 

0 

1 

5 

3 

8 

1 

1 

0 

5 

6 

9 

2 

1 

0 

5 

4 

9 

2 

6 

5 

0 

2 

4 

4 

= 

5 

5 

0 

2 

4 

4 

3 

6 

2 

0 

5 

2 

3 

4 

2 

0 

5 

2 

10 

9 

4 

5 

0 

7 

8 

9 

4 

5 

0 

7 

1 

2 

4 

2 

7 

0 

1 

2 

4 

2 

7 

0 

(2) 

0 

1 

5 

3 

8 

1 

0 

1 

5 

3 

8 

1 

1 

0 

5 

4 

9 

2 

1 

0 

5 

4 

9 

2 

5 

5 

0 

2 

4 

4 

* 

5 

5 

0 

2 

4 

4 

3 

4 

2 

0 

5 

2 

3 

4 

2 

0 

5 

2 

8 

9 

4 

5 

0 

7 

8 

9 

4 

5 

0 

7 

1 

2 

4 

2 

7 

0 

1 

2 

4 

2 

7 

0 

k<2) 


.»<4> 


a<8>-  a<5> 


Suppose  we  wish  to  find  a  minimum  cost  path  from  v^  to  v^.  The  1,5 
entry  of  A''*'  is  8  so  there  exists  a  path  of  cost  8  and  this  cost  is 
minimum.  Now  it  is  clear  that  if  there  is  a  vertex  v  such  that  the  l,j 

entry  plus  the  j,5  entry  is  8,  that  v  must  lie  on  some  minimum  cost  path. 
Consider  v,  for  instance.  The  1,4  entry  is  3  and  the  4,5  entry  is  5. 
Thus,  lies  on  a  minimum  path.  That  portion  of  a  minimum  path  connect¬ 

ing  v^  to  v-  which  is  bounded  by  v^  and  v^  must  itself  be  minimum  cost. 
Similarly  for  the  portion  between  v,  and  v  .  By  repeating  the  above  steps 
we  find  that  v^v^v^v^  is  a  minimum  dost  path  from  v^  to  v  and,  moreover, 
it  is  the  only  one.  In  general,  though,  minimum  cost  patns  are  not  unique. 
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If  we  assign  a  unit  cost  to  each  branch  of  a  network,  then  the  algorithm 
just  described  will  yield  minimum  length  paths.  This  method  differs 
from  other  algorithms  in  that  the  two  points  between  which  a  path  is  to 
be  found  need  not  be  selected  in  advance, 
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APPENDIX 


This  appendix  contains  detailed  proofs  of  all  the  theorems  in  the  text. 

Theorem  1 :  The  boundary  of  a  path  b^,,..bn  from  node  A  to  node  B  is  the 

set  of  nodes  A  and  B. 

Proof :  Let  us  calculate  the  boundary  by  Adding  the  boundaries  of  each  branch 
b^.  Notice  that  b  and  b  . ( i— 1 , 2, . . . n- 1)  are  coincident  and  hence  these 
intermediate  nodes  cancel  out  by  the  modulo  2  addition,  leaving  nodes  A  and 
B  as  the  sum. 

Theorem  2:  The  matrix  H'  (an  n  x  m  matrix)  can  be  brought  into  semi- 

diagonalized  form  by  using  the  following  two  operations: 

I.  The  addition  of  one  row  to  another 
II.  The  interchange  of  two  columns 

Proof :  In  the  first  row  there  is  at  least  one  non-zero  entry  due  to  the  in¬ 

dependence  of  the  rows.  This  can  be  brought  to  the  11  position  by  II.  Then 
the  entries  21,31,  ...  nl  can  be  made  zero  by  the  appropriate  addition  of  the 
first  row  to  the  second,  third,  ...  nth  rows.  Continuing  this  with  the  second 
row  and  so  on,  all  non-zero  entries  below  the  diagonal  are  made  zero.  Then  starting 
with  the  last  row  reverse  the  process  and  eliminate  all  non-zero  entries  above  the 
diagonal , 

Theorem  3:  If  H*  is  the  semi-diagonalized  form  of  H*  and  x*  is  obtained 

from  x  by  the  corresponding  operations  of  II  and  n*  is  obtained  from  n'  by 
the  corresponding  operations  of  I,  then  the  solutions  to  the  matrix  equation 

H*  .  x*  =  n* 

are  equivalent  to  those  obtained  from  Equation  7. 

Proof :  Let  us  express  the  matrix  Equation  7  in  the  form  of  a  system  of  simul¬ 
taneous  equations  and  analyze  the  effects  of  operations  I  and  II  on  this 
system: 


11*1 


+  h' 


,  x 
lm  m 


n 


1 


h'  x  +  . , 
nl  1 


+  h'  x  =  n ' 
nm  m  n 


The  uses  of  II  transform  this  system  to 


111 


x,  +. , .+  h’ 


l(i+l)  i+1  li  i 


+h*  x.+. .  .4h '  z  -n' 


lm  m 


h1  x  +. 
nl  1 


•4hn(i+l)Xi+l4h' 


. x.+. . ,4h *  x  «n ' 
ni  i  nm  m  n 
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The  above  is  an  instance  of  one  application  of  II,  but  in  general  the 
solutions  remain  the  same*  Let  us  assume  an  instance  of  one  application 
of  I,  where  the  ith  row  has  been  added  to  the  jth  row.  Then  all  rows 
remain  the  same  except  for  the  jth  row  which  becomes: 

(h,1l'rtlj,'l)xl+  •••  +<h’i 

Let  us  show  that  the  solutions  to  the  old  system  and  those  of  the  new  sys¬ 
tem  are  identical.  Assume  (b^,...b^)  to  be  a  solution  of  the  old  system, 
then  in  the  new  system  all  rows  other  than  the  jth  are  obviously  satis¬ 
fied.  The  jth  row  becomes 


+(h'  *'  )b  - 

im  Jm  m 


•••  h,l.b„>*<h'1lbl+  hV„>- 


■*1*— 


due  to  the  solution  satisfying  the  old  system,  and  hence  it  satisfies  the 
new  system.  Now  assume  that  (b^,  ...  b  )  satisfies  the  new  system,  then 
all  rows  of  the  old  system  other  than  tfie  jth  are  obviously  satisfied. 

The  ith  row  and  jth  row  of  the  new  system  are 


h'iibi+ 


+h'  b 
im  m 


(h’  14h1 41)b.+  ...  +(h'  +h*  )b  «  n’+n* 

il  jl  1  N  im  jm  m  i  j 


Take  the  sum  modulo  2  of  these  two  equations  and  we  get 


and  hence  the 
the  proof. 

Theorem  4:  If 


h*  b.+  ...  4h’  b  *  n* 
jl  1  jm  m  j 


Jth  row  of  the  old  system  is  also  satisfied. 


b  is  a  set  of  branches,  then  the  equation 
H  •  b  *  0 


This  completes 


is  satisfied  if  and  only  if  b  is  a  loop  or  a  set  of  loops. 

Proof :  If  b  is  a  loop  or  a  set  of  loops,  then  in  each  loop  the  boundary 
operation  determines  each  node  (that  successive  branches  are  coincident  at) 
twice,  hence  all  these  nodes  cancel  modulo  2  and  the  boundary  is  zero. 
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Assume  that  a  set  of  branches  has  a  zero  boundary.  Choose  one  branch 
and  at  one  end  there  must  be  at  least  one  other  branch  coincident  or  the 
boundary  could  not  be  zero.  The  same  holds  for  the  second  branch  and  so 
on,  since  this  cannot  go  indefinitely  we  must  finally  return  to  the  first 
branch,  hence  we  have  a  loop.  If  this  does  not  exhaust  all  the  branches 
then  repeat  the  process  with  one  of  the  remaining  branches  and  at  most 
we  have  a  set  of  loops. 

Theorem  5 


If  b  is  a  set  of  branches  and  n  is  a  set  of  two  distinct  nodes 
A  and  B  then 

H  *  b  =  n 

if  and  only  if  either  b  is  a  path  from  A  to  B  or  b  is  the  sum  of  a 
path  k  from  A  to  B  and  loops  where  k  cb  >  ^cb, .  . .  b . 

Proof : 


If  b  is  a  path  from  A  to  B  then  H*b  =  n.  If  b  is  the  sum  of  a  path 
k  and  loops  then 


b  ■  k  +  i,  +  . .  .  +  i 
1  r 

and  H*b  -  H- (k+i.+  . , .  +i  )  =  H-k+H-i,+  ...  +H •£ 

1  r  1  r 

^n  +  o-f. *fo  ~  n 

Let  us  assume  that  H*b  =  n.  At  least  one  branch  is  incident  to  A,  at 
its  other  end  it  is  either  incident  to  B  or  coincident  to  some  other  branch 
(or  else  this  node  would  be  included  in  the  boundary).  Continuing  this  pro¬ 
cess,  some  branch  will  be  incident  to  B.  If  no  more  branches  remain  in  b 
then  we  have  a  path  from  A  to  B.  If  some  branches  remain  we  choose  one 
of  them  and  it  must  be  coincident  to  at  least  one  other  branch  at  one  of  its 
ends,  and  so  on  we  find  that  the  remainder  constitutes  a  loop  or  a  set  of 
loops. 

Theorem  6 

If  A  set  of  branches  b  which  is  a  solution  to  Equation  (6  )  contains 
loops  where  i^cb , , . . i^cb ,  then  there  exists  a  path  p  such  that 

it  is  also  a  solution  to  Equation  (6  )  and  pcb. 

Proof : 

Assume  that  a  set  of  branches  b  which  is  a  solution  to  Equation  (6  ) 
contains  loops  where  i^cb , .  . .  I^cb .  Choose 

p  a  b  +  £  +  .  .  .  I 

In 
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Then 


pcb  since  i,c^,...icb  and  also 
i  n 


H* P  -  H-(b+j  +  ...  +/n)-H‘b+H'i1+..,4«.in 
-  H-b+o+  ...  +o  -  H-b  «  n 

hence  p  is  a  path. 

Lemma  1: 

If  A,  fi,  G  are  matrices  whose  entries  are  real  numbers  or  the  indeter¬ 
minate  *  then  A*(B*C)  -  (A*B)*C. 

Proof : 


Let  A  * 

f‘u) 

i,J 

-  1  ... 

n 

B  * 

M 

Ui 

-  1  ... 

n 

C  - 

M 

i»j 

*  1  .  .  . 

n 

M  -  B*C  - 

r4 

i,j 

-  1  ... 

n 

L  -  (A*B ) ■ 

(hi) 

i.J 

-  1  .  .  . 

.  n 

Then  A*(B*C)  - 

A*M  - 

■iji 

*  N 

•f'uj 

Where  r . .  *  Min  1 

k  t 

raik%y 

l- 

(A*B)xC 

-  L*C 

’hJ'M  'f;«5 

BUt  ^kj  “  ““  fak'+V  j 


,  where  r,  «  + 


>) 


ij 

,  therefore 


k  ( 


ik  kj 


J 


r  ■  Min 
lj  k 


^ik^kjj  “  |*ik^n  j^bkk,+ck’ j]J  ■ 


For  some  kQ.  r  -  aik  +  ^  and  for  aome  H'0.  Hk  j  -  \  k,  +  \ 

J  O  qj  o  o  o  o 


Therefore,  ry  »  *  +  ^  +  V  J 

o  o 


(*lk  +  bk  k'  J  +  ck'  j. 
o  o  o  o 


However,  (a^  +  bfc  Rl  )>  X  ,  and  Alk,  +  ck,  .  >r^  therefore 

0  0  o  o  oJ 


ik  '  1  k'  '  -  ~ ik •  -  ~ik'  '  k’  j  -  ij 

o  o  o  o  o  o 

r^  >  r^.  By  reversing  the  proof  we  can  show  r^>r^  hence  r 


iJ  -  ij 

which  proves  the  theorem. 


ij  ij 


Lemma  2: 


A<W  .  *«>.  *<"«)• 


We 


first  show  that  A*A^^«  A^^*A  *  ^(^+1)  For  k  *  1,  it  is  true. 
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then  A*(A(k'+1)-A*(Ak'*A)  by 


c  c  ,  , .  A*.(k')  a(u')*a  a(u'+1) 

Suppose  for  k=k\  A*A  =A  *A=A  \  „  , 

definition.  By  Lemma  1,  A*(A^C  *A)*(A*A^  )*k=k^  +^*A.  This  induction 


proves  the  first  assertion.  Now  Av  „  „  ,  ,  __ 

tinuing  this  rearrangement  we  get  A^C)*A ^^«A , 


Theorem  7. 


Q.E.D. 


Let  A  be  the  matrix  associated  with  the  network  N  and  the  cost 


function  f. 


Let  a_  ^  be  the  ij  entry  of  A^. 


Then  a , . 

iJ 


00 


is  the 


least  cost  entailed  by  any  path  consisting  of  <k  branches  joining  v.  to 

Y 

Proof: 

Let  X  <h)  be  the  least  cost  entailed  by  any  path  consisting  of  <h 
ij  {i\ 

branches  and  joining  v,  to  v.,  It  is  clear  that  X..  *  a,.,  where  a,, 

i  j  ij  ij  ij 

is  the  i,J  entry  of  A,  for  there  is  but  one  path  joining  to  v 

having  <1  branch.  Either  it  is  the  path  having  no  branches,  in  whicn  case 
«a.  ■  oo  or  it  is  the  path  consisting  of  the  branch  (v  ,v  )  in  which 
case  X^j  w  Let  A  be  the  matrix  associated  with  the  network  N 

and  the  cost  function  f.  Let  a^j^  be  the  ij  entry  of  A^tnake  induction 
assumption,  etc..  Let  X. . ^  be  the  cost  entailed  by  the  cheapest  path  from 
to  vi  having  <k  branches.  Let  p  be  such  a  path.  Either  p  consists 

--  j -  —  _ranch  (v^,v  )  or  there  is  a  point  v^  on  p  intermediate  to  v^ 

and  Min^a^^k  ^  assumption  k>l.  Hence  for  p*j 

a  is  the  cost  of  a  cheapest  path  from  v,  to  v  with  <k-l  branches  .* 

fli*J  (k)M  ^  (k)  \t _ _  r- _  _  1 1  _r  _ _ _  i_  _  j _ i_  1.1 _ i_  i  _  i _ _ 

ij  (h)  ,  (h)  .  .  .  ,  .  A(h) 

where  a^j  1S  t^e  L>j  entry  in  A 


\  j ^ow  suPPose  f°r  flii  integers  h  such  that  l<h<k,  we  have 

X.  .(h)  *  where  a  is  the  i,j  entry  in  A^.  Consider  X. 

ij  ij  ij  ,J  7  ij 

There  is  a  path  p  connecting  v^  to  v.  and  having  <k  branches  such  that 

X.  .(U>  is  the  cost  of  p.  We  assume  p  ilas  at  least  one  branch.  There  are 
ij  * 

two  cases  to  consider: 


Case  I: 

p  consists  of  the  branch  joining  v.  to  v  .  Then  X, .  X  ^ 

,  (k-1),  ,  (k-1)  (k-1)  1  /  (k-1) 

x.  .  '  +  o  ■  X.  .  '+  a. .  ■  a. .  +  a,  ,>Min  (a.,  +a, 

1J  1J  4m  ,J|r  lh  ij 

If  for  some  h  £  j ,  we  have  X. . ^K^>a ,u i^+a,  , 

ij  in  nj 


a. 


) 


then  we  may  conclude  there 


is  a  path  of  <k-l  links  going  from  v.  to 


and  then  going  to  v  by 
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one  more  branch  whose  cost  is  less  than 
mallty  of  the  cost  of  p. 

Case  2: 


vij 


(10 


This  contradicts  the  mini- 


There  is  a  third  point  v .  intermediate  Co  and  v.  through 

which  p  passes.  Let  t  be  the  number  of  branches  on  thatJpart  of  p 

(10  _  ^  (t)+x(k-t)_ 


connecting 
a 


V  j  to  V..  iU«=LI  Ww  W..WW  /V,  .  -  /v.. 

.  (fc)+a  0c“t)>A .  S**  •  Again  as  on  Case  1,  if  X  H)>a,  .  ^ 

SL. 


Then  we  see  that  X. 


there  exists 


Q.E.D. 


Theorem  0: 


If 

for  all 
Proof: 

The 


k>k' 

i  >j  • 


then  the  i,j  entry  of  A^<  the  i,j  entry  of  A^  ^ 


i,j  entry  of 


A<k’>  is 


(k') 


consider  longer  paths  (consisting  of  more  branches)  connecting  vJt  v1(  we 


may  come  up  with  a  cheaper  path,  in  which  case 
not,  in  which  case  X^^®  j  ^  ^  * 


by  the  previous  theorem. 

tine 

(k’? 


If  we 


x  00<x 

ij 


±>  y 

or  we  may 


Q.E.D. 

In  view  of  lemma  2  and  theorem  2,  if  we  are  dealing  with  a  network  of 
nodes,  the  most  matrix  compositions  we  need  to  do,  in  order  to  determine 


all  minimum  costs  between  nodes,  is  logjn+l. 

(2)  A (4) 


is  to  perform  the  "doublings"  A 


The  procedu 

<2l5g 2nVt> 

, . . .  .A 


re  for  getting  A 


(ti-1) 


An^(n’l). 


>-l) 


necting  v . ,v  . 

a  <n-V  j 
ij  * 

for  all 


we  can  easily  determine  a  minimum  cost  path  con- 
oo  there  are  no  paths  connecting 


Once  we  have  A 

If  A..'  _r*  oo  there  are  no  paths  connecting  v,v 

then  look  for  a  node  v^  such  that  h^j,  h^/^i  a^ 

1  (n-l) .  (n-1)  _  (n-1) 

a  i  u  "t*  au  4 


h/»j  and  for  which  a^ 


.  hi^ 

node  exists  then  the  branch  connecting  v  to  v  is  already  a  minimum  cost 
does  exist  then  there  is  a  branchJof  the  network  N  connect  - 

In  that  case  we  repeat  the  procedure  to  find  a  minimum  cost 


i  j 


V  v  If 


If  no  such 


path.  If 
ing  v.  to  v 


j* 

path  connecting  v^  to  vh^»  etc* 

Illustrations: 

Page 


9 

- 

Dwg. 

IB-10,027 

10 

- 

11 

U- 10, 026 

19 

- 

It 

ID-10,028 

20 

- 

11 

ID-10,029 

21 

- 

II 

ID- 10, 030 

22 

- 

11 

ID-10,031 

3.  On. 

D.  M.  Liss 

- - 


S.  Olcada 


1-6 


ffC.  J.  Kolker- 


TM-3421 


REFERENCES 


1.  B-TOG:  C.  Berge  :MTheory  of  Graphs,*1  J.  Wiley  &  Sons,  New  York,  1962. 

2.  L-apa:  C.Y.  Lee:  “An  Algorithm  for  Path  Connections  and  Its  Appli¬ 
cations11,  IRE  Trans.  PG. EC-10. 3  September  1961,  346-65. 

3.  M-SPTM:  E.  F.  Mcore:  “The  Shortest  Path  Through  a  Maze",  Proc. 

Inter.  Symp.  on  The  Theory  of  Switching" , Harvard  Univ,  Press,  1957. 

4.  O-tasc:  S.  Okada:  “Topology  Applied  to  Switching  Circuits",  Proc. 
Symp.  on  Information  Networks,  P.I.B.,  Brooklyn,  pp.  267-290,  April, 
1954. 

5.  V-AS:  0.  Veblen:  “Analysis  Situs",  Amer.  Math.  Soc.  Cell.  Pub., 

Vol.  V,  Part  2,  Amer.  Math.  Soc.,  N.Y. ,  1931. 


II-l 


£0  CO  - 
ft  ft  O  ft 


P  W  •P  TJ 
O  CO  to  0) 


ft  CO  ft  ft 

co  X 
f>>  *N  O  0 
d  ft  o 


d 

)  I  iu  • 

P  O  Pi 

O  CD 

Ci  T)  -P 

•p  g  p 

CO  ft 

P  R 

p  ft  o 

5  S3  ° 

>  CO 


0  l 

>  ft  ft 
•H  p  0 
bD  S 
*H  d 
W  (0  O 
•H  d  ft 
_  P  d 


ft  ft  £  D<  Cfl 


£  ft  rH  o 


CQ  CO  • 

♦h  p  a  p 


p  to  ft  ft 

O  CO  CQ  CD 


ft  CO  p  p 

CO  X 

>?  _?  o  <u 

d  ft  O 

CO  0)  0 

_  d  P  ft 

ft  ft  CO 


CO 

d 

JjOICO* 

a)  d  ft  a  d 

ft  ft  O  0 

ft  ft  0  ft  ft 

d  p  d  d 
ft  ft  CO  ft 

o  ft  p  W 

ft  ft  O 


ft  o  d  ft  p 

d  P  0  I 
O  ft  >  ft  ft 

O  0  ft  P  0 

0  H  tD  bO  0 


•H  d 

CO  co  O 

ft  d  ft 
.  ft  d 
P  CO  0 
ft  ft 


O  ft  d 

0ft 

“ 

ft  0  0 
_  ft  d 


co  d  e  0  ft 

p  0  ft  p  ft 

P  ft  d  w 
0  0  ft  ft  0 

ft  ft  g  d  0 


K  ,  O 

ftp 

d  Eh  o 


